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HELDMG : Heliothis  (Lepidoptera: 

Noctuidae)  Spp.  Damage  Subroutine  for  the 
Cotton  Model  GOSSYM 

W.M.  Thomas 


HELDMG  is  a computer  model  describing  the 
within-plant  distribution  of  damage 
caused  by  Heliothis  spp.  larvae  in 
cotton.  It  will  be  added  to  the  cotton 
model  GOSSYM  (Baker  et  al.  1983)  as  a 
subroutine  and  provide  a means  of 
integrating  a Heliothis  spp.  model  with 
GOSSYM.  This  integrated  model  can  be 
used  to  study  the  effects  of  Heliothis 
spp.  damage  on  cotton  growth  and  yield. 

HELDMG  consists  of  a command  procedure, 
INMENU,  and  five  subroutines — RDDMG, 
PREDMG,  DAMAGE,  REDIST,  and  FRMTRX. 
However,  when  a Heliothis  spp.  model  is 
integrated  with  GOSSYM,  the  command 
procedure  and  the  subroutine  RDDMG  are 
not  needed.  The  program  is  written  in 
FORTRAN  and  has  been  run  on  a VAX 
11/750.  A general  flow  chart  of  HELDMG 
is  shown  in  figure  1. 


COMMAND  PROCEDURE  INMENU 


A command  procedure  file  is  included  to 
illustrate  how  the  model  uses  different 
types  of  input  data.  It  is  executed  at 
the  beginning  of  the  simulation  and 
prompts  for  Heliothis  damage  data.  A 
choice  is  required  from  the  following 
menu: 

Input  Heliothis  data 

1.  Heliothis  model 

2.  Scouting  report 

3.  No  damage 


Research  entomologist.  Cotton  Production 
Research  Unit,  Agricultural  Research 
Service,  U.S.  Department  of  Agriculture, 
P.0.  Box  2131,  Florence,  SC  29503. 

This  study  was  conducted  by  the 
Agricultural  Research  Service,  USDA,  in 
cooperation  with  the  South  Carolina 
Agricultural  Experiment  Station. 


Figure  1 

General  flow  chart  of  HELDMG. 
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Option  1:  If  you  use  information  from  a 
Heliothis  model,  another  menu  will  be 
presented.  From  this  second  menu  you  can 
either  (A)  enter  the  number  of  larvae  per 
acre  or  (B)  enter  the  number  of  fruit 
damaged  by  small  and  large  larvae.  These 
data  are  obtained  from  the  Heliothis  spp. 
model  MOTHZV  (Hartstack  and  Witz  1983). 

If  you  choose  option  A,  then  the  day  of 
the  year  of  infestation,  physiological 
age  of  the  larva,  and  the  number  of 
larvae  per  acre  must  be  entered.  If  you 
choose  option  B,  then  the  day  of  the  year 
of  infestation  and  number  of  fruit 
damaged  per  acre  by  small  larvae  and  by 
large  larvae  must  be  entered. 


classes — small,  medium,  and  large 
squares,  flowers,  and  small,  medium,  and 
large  bolls. 

The  number  of  damaged  fruit  is 
distributed  among  these  age  classes, 
depending  on  the  type  of  input  data.  If 
the  data  are  from  a Heliothis  model,  the 
damage  is  distributed  based  on  preference 
and  stratification  coefficients  derived 
by  Wilson  and  Gutierrez  (1980) . These 
coefficients  describe  the  changing 
preferences  of  Heliothis  larvae  for 
different  aged  fruits  as  the  age 
distribution  of  the  larvae  and  fruit 
within  the  cotton  plant  changes. 


Option  2:  Data  from  scouting  reports 
concerning  Heliothis  damage  can  also  be 
entered  into  the  model.  These  data  are 
day  of  year  of  infestation,  percent 
damaged  squares,  and  percent  damaged 
bolls. 

Option  3:  This  option  runs  GOSSYM  with  no 
Heliothis  damage. 

A temporary  work  file  ( INPUT. DMG) 
containing  the  Heliothis  damage 
information  is  created.  Data  in  tnis 
file  are  used  by  the  subroutines. 


SUBROUTINES 


RBDMG 

For  each  simulated  day,  INPUT. DMG  is 
read.  If  the  number  of  larvae  per  acre 
has  been  entered,  and  if  damage  is  to 
occur  today,  then  the  number  of  fruit 
damaged  today  is  calculated  using  an 
equation  derived  by  Townsend  (1973). 

HZBUG  = 0.035  + (0.093  * BUGAGE) 

where: 

HZBUG  = number  of  fruit  damaged  today 
BUGAGE  = physiological  age  of  larva 

PREDMG 

The  cotton  plant  is  inventoried,  and  the 
fruiting  forms  are  grouped  into  seven  age 


^ij  AiSi  jCi/l  (AfSi  jCf  j ) 
where: 


ij 


proportion  of  attacks  b 
instar  larvae  against  i 
fruit  age  class  (CKP<1) 


A^  = proportion  of  total  fruit  in 
fruit  population  belonging  to 
i^  age  class  (0<A<1) 


S-j^j  = relative  preference 

(0<S<1)  that  instar 
larvae  have  for  i1-*1  fruit  age 
class 


C-j^j  = relative  within-plant 

stratification  coefficient 
(0<C<1)  of  jt-k  instar 
larvae  for  i*-*1  fruit  age  class 


The  number  damaged  in  each  age  class  is 
totaled  for  small  larvae  (instars  1-3) 
and  large  larvae  (instars  4-5) . This 
number  is  passed  to  the  DAMAGE  subroutine 
and  is  distributed  within  the  cotton 
plant . 


If  percent  square  and  boll  damage  are 
entered,  then  the  number  of  fruit  on  the 
plant  is  estimated  by  GOSSYM.  The 
percent  damage  is  converted  to  number 
damaged.  This  damage  is  distributed 
among  the  age  classes  of  fruit  by 
assuming  that  the  amount  of  damage 
occurring  in  each  age  class  is  directly 
proportional  to  the  number  present.  The 


2 


total  number  of  fruit  damaged  within  each 
age  class  is  calculated  and  this  number 
is  passed  to  the  subroutine  DAMAGE. 

FRUTDM3(I)  = PR OP (I)  * NUM 

where: 

FRUTDM3(I)  = number  of  fruit  in  each 
age  class  damaged  by 
Heliothis  larvae 

PROP(I)  = proportion  of  fruit  in 
each  age  class 

NUM  = number  of  fruit  damaged 

per  acre  by  Heliothis 

DAMAGE 

Most  of  the  Heliothis  damaged  fruit 
occurs  in  the  upper  part  of  the  plant 
(Gonzalez  et  al.  1967,  Fye  1972,  Ramahlo 
et  al.  1984)  and  close  to  the  mainstem, 
with  over  80  percent  within  three 
fruiting  sites  from  the  mainstem 
(Nicholson  1975,  Ramahlo  et  al.  1984). 


Therefore,  this  subroutine  weights  the 
distribution  of  damaged  fruit  vertically 
and  horizontally  on  the  plant. 

It  is  assumed  that  the  number  of  fruit 
damaged  on  the  mainstem  and  vegetative 
stems  is  proportional  to  the  number  of 
fruit  on  these  stems.  Damaged  fruit  is 
distributed  vertically  on  the  mainstem 
based  on  the  larval  instar.  The  mainstem 
is  divided  into  three  equal  regions  and 
the  vegetative  stems  are  grouped  in 
region  4 (fig.  2).  The  amount  of  damage 
occurring  in  each  region  on  the  mainstem 
is  based  on  information  from  Ramahlo  et 
al.  (1984)  (table  1).  These  percentages 
are  used  in  the  model  to  distribute  the 
Heliothis  damage  vertically  on  the  cotton 
plant . 


— x 


X — X— 


MAINSTEM 


VEGETATIVE  STEMS 


Figure  2 

Diagram  of  simulated  cotton  plant  showing 
regions  used  in  distributing  Heliothis 
damaged  fruit.  (x  = fruiting  site) 
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Table  1 

Percent  of  fruit  damage  caused 
by  Heliothis  larvae  occurring 
vertically  per  region  on  cotton 
plant ' s raainstem 


Region  of 
mainstem 

Small 

larvae^/ 

Large 
larvae 2/ 

Scouting 

report 

Upper  1/3 

50 

39 

45 

Middle  1/3 

37 

43 

40 

Lower  1/ 3 

13 

18 

15 

After  the  number  of  damaged  fruit  is 
distributed  on  the  mainstem,  the  fruiting 
sites  on  the  vegetative  stems  are 
evaluated.  Since  a small  proportion  of 
the  total  number  of  fruit  is  on  the 
vegetative  stems,  no  attempt  is  made  to 
weight  the  distribution  of  the  damaged 
fruit.  The  number  of  fruit  to  be  damaged 
in  each  age  class  is  compared  with  the 
number  present.  If  there  is  more  fruit 
to  be  damaged  than  available,  then  all 
fruit  is  damaged  and  the  remainder  is 
stored.  If  enough  fruit  is  available, 
then  the  percent  of  each  fruit  to  be 
removed  is  calculated. 


1/ Instars  1-3. 

2/lnstars  4-5. 

Once  the  number  of  fruit  to  be  damaged 
within  each  region  on  the  mainstem  is 
calculated,  the  damage  is  distributed  on 
each  fruiting  branch  based  on  the 
distance  of  the  fruiting  site  from  the 
mainstem.  Ramahlo  et  al.  (1984)  found 
that  approximately  54  percent  of  the 
damage  occurred  at  the  first  site  on  the 
fruiting  branch,  23  percent  at  the 
second,  12  percent  at  the  third,  5.5 
percent  at  the  fourth,  and  5.5  percent  at 
the  fifth.  The  amount  of  fruit  damage  at 
each  position  on  the  branch  is  calculated 
by— 


VDMG(IC) 
where : 


s CREG ( 4 , I C ) / VEG ( I C ) 


VDMG(IC)  - percent  of  fruit  to  be 
damaged  in  each  age 
class  on  vegetative 
stems  because  of 
Heliothis  larvae 

CREG(4,IC)  = number  of  fruit  to  be 
damaged  on  vegetative 
stems  for  each  age  class 

VEG(IC)  = number  of  fruit  in  each 
age  class  on  vegetative 
stems 


POSDMG(IR,M,IC) 


where : 

POSDMG(IR,M,IC) 


CREG(IR,IC) 


BPER(M) 


CREG(IR,IC)  * If  any  fruit  remains  to  be  damaged,  it  is 

BPER(M)  equally  distributed  among  all  fruiting 

sites . 


ADDMG 

number  of  fruit  per 
region  per  position  where: 
per  age  class 

damaged  by  Heliothis  ADDMG 


number  of  fruit 

damaged  per  region 

per  age  class  REMAIN 


amount  of  damage 
occurring  at  each 
fruiting  site  (54% 
position  1,  23% 
position  2,  12% 
position  3,  and 
5.5%  positions  4 
and  5) 


S FRUIT 

ACTDMG 


REMAIN/ (SFRUIT  - ACTDMG) 


amount  of  additional  fruit 
to  be  damaged  on  each 
fruiting  site 

number  of  fruit  remaining  to 
be  damaged 

number  of  fruit  susceptible 
to  Heliothis  damage 

actual  number  of  fruit 
damaged 
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The  percent  removed  from  each  fruiting 
site  on  the  mainstem  is  calculated  by — 


HZDMG 


where: 


HZ  DM3 


= (POSDMG(IR,M,IC)/TCLS 
(IR,M,IC) ) + ADDMG 


percent  of  fruit  to 
be  removed  from 
fruiting  site  due 
to  Heliothis  feeding 


POSDMG( IR,M ,IC)  = 


TCLS( IR,M,IC) 


number  of  fruit  per 
region  per  position 
per  age  class  to  be 
damaged  by  Heliothis 

number  of  fruit  per 
region  per  position 
per  age  class 
susceptible  to 
damage 


ADDMG  = amount  of 

additional  fruit  to 
be  damaged  on  each 
fruiting  site 

HZDMG  is  used  to  decrement  the  amount  of 
square  and  boll  weight  and  fraction  of 
each  fruiting  site. 

REDIST 

This  subroutine  redistributes  any 
remaining  damaged  fruit  within  a region. 

FRMTRX 

This  subroutine  is  called  at  the  end  of 
the  simulation  and  prints  the  percent 
damage  occurring  at  each  fruiting  site. 


AGE(3,30,5) 

AGEB0L(3 ,30,5) 

BEGIN (4) 

BOLABZ 

B0LWGT(3 ,30 , 5) 
BPER(5) 

BUGAGE 

BUGDM3 (7) 

BUGGY 

BURRN 
CD  M3 

CHZDM3(3,30,5) 

CLASS (3,30,5) 


DEFINITION  OF  TERMS 


COTXX 


ACTDM3 


ADDMG 


actual  number  of  fruit 
damaged  per  plant 

amount  of  additional 
fruit  to  be  damaged  on 
each  fruiting  site 


CREG(4 ,7) 


DENOM 


DM3  (7) 


= age  of  each  node  (days 
after  node  initiation) 

= age  of  boll  (days  after 
blooming ) 

= first  node  of  region 

= number  of  boll 
abscissions 

= boll  weight  (grams) 

= amount  of  damage 
occurring  at  each 
fruiting  site 

= physiological  age  of 
larva 

= number  of  fruit  in  each 
age  class  on  mainstem 
damaged  by  Heliothis 

= a logical  constant  used 
to  signal  if  Heliothis 
damage  has  occurred  today 

= burr  nitrogen 

= counts  number  of  damaged 
fruit 

= cumulative  loss  of  fruit 
at  node  due  to  Heliothis 
damage 

= age  class  of  each  fruit 
(1  = small  square,  2 = 
medium  square,  3 = large 
square,  4 = flower,  5 = 
small  green  boll,  6 = 
medium  green  boll,  7 = 
large  green  boll) 

= weight  in  grams  of  open 
bolls 

= number  of  fruit  damaged 
per  region  per  age  class 

= denominator  used  in 
calculating  probability 
of  larval  attack 

= number  of  fruit  damaged 
per  age  class 
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DT0TAL(7 ) 

number  of  fruit  in  each 
age  class  on  mainstem 

N 

counter 

NBRCH 

= 

number  of  fruiting 

END (4) 

= 

last  node  of  region 

branches 

FCODE( 3,30 ,5) 

= 

fruit  code  (1  = square, 

NFBR(l) 

= 

number  of  fruiting 

2 = green  boll,  3 = open 
boll,  4 = abscised,  3 = 

branches  on  mainstem 

square  marked  for 

NNID 

= 

node  number  on  fruiting 

abscission,  6 = boll 

branch 

marked  for  abscission) 

NNOD(3,30) 

= 

number  of  nodes  on 

FFRUT(3,30 ,5) 

= 

fraction  of  fruit 

fruiting  branch 

remaining  at  node 

NSTR 

= 

larval  instar 

FRUTDMG ( 7 ) 

= 

number  of  fruit  in  each 
age  class  damaged  by 

NUM 

= 

number  of  larvae  per  acre 

Heliothis  larvae 

NVBRCH 

= 

number  of  vegetative 

FRUTLOS 

= 

fraction  of  fruit  lost 

branches 

GBOLWT 

= 

total  green  boll  weight 

PD MG (5, 7) 

= 

probability  of  damage 
occurring  in  each  fruit 

GBP 

= 

percent  of  green  boll 

age  class  by  each  larval 

damage 

instar 

HZBUG 

= 

number  of  fruit  damaged 

PLTDIC 

= 

number  of  fruit  damaged 

by  Heliothis 

per  plant 

HZDMG 

= 

percent  of  fruit  to  be 
removed  from  fruiting 

POPPLT 

= 

plant  population 

site 

P0SDM3(4,5,7) 

number  of  fruit  per 
region  per  position  per 

IBUGDAY 

= 

day  of  year  damage  is  to 

age  class  damaged  by 

occur 

Heliothis 

IC 

= 

fruit  age  class 

PQFLR 

= 

weight  lost  from  fruit 
due  to  petal  shed  after 

ICLS 

= 

fruit  age  class 

blooming 

IOPTION 

= 

menu  selection  relating 

PREF(5 , 7 ) 

= 

larval  food  preference 

to  type  of  input  data 

values  for  various  age 
fruit 

IPOS  (4) 

' 

maximum  number  of 
positions  on  fruiting 

PREFA(2 , 7 ) 

= 

larval  food  preference 

branch 

values  for  various  age 
fruit  for  small  (instars 

IR 

region 

1-3)  and  large  (instars 
4-5)  larvae 

LGLARV ( 7 ) 

number  of  fruit  in  each 
age  class  damaged  by 

PR0P(7) 

= 

proportion  of  fruit  in 

large  larvae  (instars 

each  age  class 

4-5) 

PVEG 

= 

proportion  of  fruit  in 

LOWER 

lower  1/ 3 of  mainstem 

each  age  class  on 
vegetative  stems 

MID 

= 

middle  1/ 3 of  mainstem 
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RDMG 

REG (4,7) 
REGDM3  (4,7) 

REMAIN 

REMN 


RPOS 

SEEDN 
S FRUIT 

SMLARV(7) 


SQP 

SQRWT( 3,30 ,5) 
SQRZ 

STRAT( 5 , 7 ) 

STRATA (2, 7) 


TCLASS (7) 
TCLS  (4,5,7) 


TOTGB 


- number  of  fruit  TOTSQR 

remaining  to  be  damaged 

UPPER 

= number  of  fruit  per  age 

class  per  region  VDMG(7) 

= number  of  fruit  damaged 
per  region  per  age  class 


= number  of  squares 

- upper  1/ 3 of  mainstem 

= percent  of  fruit  in  each 
age  class  on  vegetative 
stems  damaged  by 
Heliothis  larvae 


- number  of  fruit  VEG(7) 

remaining  to  be  damaged 


= number  of  fruit  left 
over  to  be  damaged  that 
will  be  added  to  next 
position  on  fruiting 
branch 

= number  of  fruiting  sites 
on  fruiting  branch 


VEGP 


WTLOS 

XDIFF 


- seed  nitrogen 


= number  of  fruit  in  each 
age  class  on  vegetative 
stems 

= proportion  of  fruit  in 
each  age  class  on 
mainstem 

= weight  loss 

= number  of  fruit 

remaining  to  be  damaged 


= number  of  fruit  EXAMPLE  OF  INPUT 

susceptible  to  Heliothis 
damage 


Scouting  report 

data 

= number  of  fruit  in  each 

age  class  damaged  by 
small  larvae  (instars 

Day  of  year  % 

square  damage 

% boll  damage 

1-3) 

169 

4.0 

0 

17  6 

3.0 

0 

= percent  of  square  damage 

181 

1.0 

0 

188 

1.0 

0 

= weight  of  squares 

196 

4.0 

0 

201 

4.2 

4.2 

= total  squares 

209 

11.0 

5.3 

216 

17.0 

11.0 

- age  specific  larval 

222 

7.0 

6.0 

stratification 

228 

2.0 

6.0 

coefficients 


= age  specific  larval 
stratification 
coefficients  for  small 
and  large  larvae 

= number  of  fruit  for  each 
age  class 


= number  of  fruit  per 
region  per  position  per 
age  class  susceptible  to 
damage 


= number  of  green  bolls 
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EXAMPLE  OF  OUTPUT^/ 


Percent  fruit  damaged  at  each  location 


Mainstem  (Kl) 


Vegetative  branch  1 (K2) 


Ml  M2  M3  M4  M5  Ml  M2  M3  M4  M5 


L 1 
L 2 
L 3 
L 4 
L 5 
L 6 
L 7 
L 8 
L 9 
L 10 
L 11 
L 12 
L 13 
L 14 
L 15 


8 2 3 1 0 

8 3 10  0 

23  9 3 0 0 

22  10  7 0 0 

26  6 0 0 0 

25  13  0 0 0 

27  14  0 0 0 

17  0 0 0 0 

17  0 0 0 0 

17  0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 


20  0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

Vegetative  branch  2 (K3) 

Ml  M2  M3  M4  M5 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 


1/ From  HELDMG  using  scouting  report  data 
as  input.  (K  = vegetative  branch  number, 
L = fruiting  branch  number,  M = fruiting 
branch  node  number) 
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FLOW  CHARTS 


RDDMG 


Rewind  INPUT. DMG. 


Read  input  option. 


Are  input  data  from 
Heliothis  model? 


Input  day  of  year,  number  of 
larvae  per  acre,  and  physiological 
age  of  larvae. 


Will  damage  occur  today  ? 


Calculate  number  of  fruit  damaged 
per  acre.  This  equation  was 
derived  by  Townsend  (1973). 


If  number  of  fruit  damaged  is  less 
than  0,  set  it  to  0. 
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0 0 


0 0 


Assign  instars  (NSTR)  to  larvae  based 
on  physiological  age. 
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If  damage  will  occur  today,  set 
BUGGY  = TRUE. 

Call  PREDMG  subroutine  to  distribute 
total  number  of  damaged  fruit 
among  different  age  classes  of  fruit. 


Are  input  data  from  MOTHZV? 


Input  day  of  year,  number  of  damaged 
fruit  per  acre  caused  by  equivalent 
number  of  small  (DMGSML)  and  large 
(DMGLRGE)  larvae. 


Vill  damage  occur  today? 
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Have  small  larvae  caused  damage? 


Set  NSTR  = 1 and  HZBUG  = number  of  fruit 
damaged  by  small  larvae. 


Call  PREDMG  subroutine  to  distribute 
damage  among  different  age  classes  of 
fruit.  Set  BUGGY  = true. 


Have  large  larvae  caused  damage? 


Set  NSTR  = 2 and  HZBUG  = number  of  fruit 
damaged  by  large  larvae. 


Call  PREDMG  subroutine. 


If  BUGGY  = true,  call  DAMAGE  subroutine. 
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FALSE 


Are  input  data  from  scouting  report? 


0 


Input  day  of  year,  percent  square  damage, 
and  percent  green  boll  damage. 


Will  damage  occur  today  ? 


If  percent  square  damage  or  percent  green 
boll  damage  is  greater  than  0,  call  PREDMG 
subroutine. 


Call  DAMAGE  subroutine  to  distribute 
damaged  fruit  on  plant. 
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PREDMG 


Set  all  counters  to  0. 


Do  for  each  fruit  age  class  (I)  and 
each  larval  instar  (5). 


Set  counters  to  0. 


Assign  fruit  to  different  age  classes. 
Examine  all  monopodial  branches  (K), 
all  sympodial  branches  (L),  and  all 
fruiting  sites  on  each  sympodium  (M). 


Fruit  must  be  between  2 and  48  days 
old  to  be  susceptible  to  Heliothis  damage. 
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Classify  green  bolls. 


If  green  boll  is  less  than 
3 days  old,  then  consider  it  a 
bloom. 


If  green  boll  is  between 
2 and  1 0 days  old,  then  consider 
it  a small  boll. 


If  green  boll  is  between 
9 and  17  days  old,  then  consider 
it  a medium  boll. 


If  green  boll  is  greater 

than  16  days  old,  then  consider 

it  a large  boll. 
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Classify  squares. 


If  square  is  between  1 and  9 
days  old,  then  consider  it  a small 
square. 


If  square  is  between  8 and  17 
days  old,  then  consider  it  a medium 
square. 


If  square  is  older  than  16  days, 
then  consider  it  large. 


Accumulate  number  of 
susceptible  fruit  (SFRUIT) 
and  number  of  susceptible 
fruit  per  age  class  (TCLASS). 


Accumulate  number  of 
squares  (TOTSQR). 


Accumulate  number  of 
green  bolls  (TOTGB). 
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FALSE 


Input  data  from  Heliothis  model. 


IOPTION  = 1 

\ o / 


TRUE 


Do  for  each  fruit  age  class. 


If  no  fruit  is  available , 
return  to  main  program. 


Calculate  proportion  of  each  age 
class  of  fruit  (PROP). 


Do  for  each  fruit  age  class. 


Distribute  damaged  fruit 
among  age  classes  based 
on  preference  (PREF)  and 
stratification  (STRAT) 
coefficients  derived  by 
Vilson  and  Gutierrez  (1 980). 


(D 


© 


© 


For  each  fruit  age  class  : 


Calculate  number  of  fruit  of  each 
age  class  damaged  by  small  larvae 
(SMLARV)  and  large  larvae  (LGLARV). 


Input  data  from  MOTHZV. 


For  each  fruit  age  class  : 


If  no  fruit  is  available, 
return  to  main  program. 


Calculate  proportion  of  each  age  class 
of  fruit. 
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0 


0 


For  each  fruit  age  class  : 

Same  logic  as  Heliothis  model, 
except  preference  and  stratification 
coefficients  averaged  for  instar  1 -3 
and  instar  4-5. 

For  each  fruit  age  class  : 


For  each  fruit  age  class  : 


Calculate  number  of  fruit  to  be 
damaged  by  large  larvae. 


Calculate  number  of  fruit  to  be 
damaged  by  small  larvae. 
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Input  data  from  scouting 
report. 


If  percent  square  damage 
is  0,  skip  to  green  boll 
section. 


For  each  size  square  : 


Calculate  proportion  of 
squares  in  each  age  class. 


Calculate  number  of  fruit 
damaged  per  acre  (NUM). 


For  each  size  square  : 


Number  of  fruit 
damaged  in  each  age 
class  is  proportional 
to  number  present. 
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Same  logic  for 
green  bolls. 


For  each  size 
green  boll  : 


For  each  size 
green  boll  : 
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DAMAGE 


Set  ell  counters  to  0. 


For  each  fruit  age  class  (1),  each 
region  (K),  and  each  fruiting  site  (L) 
on  sympodium  : 


Set  counters  to  0. 


Divide  plant  into  thirds  based  on 
number  of  mainstem  nodes. 


Consider  each  stem. 

1 = mainstem. 

2 & 3 = vegetative  stems. 


Regions  1 -3  are  for  mainstem  only 
and  region  4 is  for  both  vegetative 
stems. 
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A 


Assign  beginning  and  ending  node 
for  each  region. 
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This  section  counts 
number  of  fruit  in  each 
age  class  available  for 
damage. 

For  each  fruiting  site  on 
sympodium  : 


Count  number  of 
fruit  per  age  class  on 
mainstem  (DTOT  AL) 
and  store  maximum 
number  of  positions  per 
region  (IPOS). 


Count  number  of 
fruit  on  vegetative 
stems  (VEG)  and  set 
maximum  number 
of  positions  to  5. 

Count  number  of 
fruit  available  to 
damage  (SFRUIT), 
number  of  fruit  per 
age  class  in  each  region 
(REG),  number  of  fruit 
per  age  class  per  region 
in  each  position  (TCLS). 
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0 


For  each  fruit  age  class  : 


Calculate  proportion  of  fruit  on 
mainstem  (VEGP)  and  vegetative  stems 
(PVEG)  for  each  age  class. 


Information  from  Heliothis  model. 
Small  larvae  - instars  1-3. 


Number  of  fruit  to  be  damaged  per  age  class  is 
converted  to  per  plant  basis  (PLTDMG)  and 
number  to  be  damaged  is  weighted  vertically 
for  each  region  (CREG)  based  on  information 
from  Ramahlo  et  al.  (1  984). 


Large  larvae  - instars  4-5. 
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Information  from  scouting  reports. 


If  no  fruit  of  this  age  class  is  to  be 
damaged,  skip  to  next  age  class. 


Number  of  fruit  to  be  damaged  per  age 
class  is  converted  to  per  plant  basis. 


Number  to  be  damaged  is  weighted  vertically 
for  each  region. 


Count  number  of  fruit  to  be  damaged 
per  age  class  on  mainstem  (BUGDMG). 
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For  each  fruit  age  class  : 


If  no  fruit  of  this  age  class  is  to  be 
damaged,  skip  to  next  age  class. 


If  number  of  fruit  in  this  age  class 
to  be  damaged  is  greater  than  or  equal 
to  number  available,  then  set  DMG 
for  that  age  class  to  1 .0  and  move  any 
remaining  fruit  to  be  damaged  to 
vegetative  stems. 


For  each  region  : 


If  no  fruit  is  in  this  age  class 
vithin  this  region,  then  move  number 
to  be  damaged  to  next  region. 


If  number  of  fruit  in  this  age  class 
to  be  damaged  is  greater  than  or  equal 
to  number  available,  then  set  REGPMG 
for  that  age  class  to  1 .0  and  move  any 
remaining  fruit  to  be  damaged  to  next 
region. 
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If  more  fruit  is  available 
than  needs  to  be  damaged,  then 
consider  each  fruiting  site. 


For  each  fruiting  site  on 
sympodium  : 


If  fruit  at  this  site  is  in  correct 
age  class  and  this  age  class  is  to  be 
damaged,  then  calculate  percent  to 
remove  from  this  site  (POSDMG).  Damage 
is  weighted  horizontally  based  on  distance 
from  mainstem  (BPER).  Weighting 
is  based  on  information  from  Ramahlo  et  al. 
(1984). 


If  number  of  fruit  in  this  age  class 
to  be  damaged  is  greater  than  or  equal 
to  number  available,  then  set  POSDMG 
equal  to  number  available  (TCLS)  and 
move  any  remaining  fruit  to  be  damaged  to 
next  position  on  fruiting  branch. 
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Count  number  of  fruit  damaged. 


If  not  enough  fruit  has  been 
damaged  from  this  region, 
then  call  RED  1ST  to  redistribute 
fruit  by  position. 


For  each  fruit  age  class  : 

Check  vegetative  stems 
for  damage. 


If  number  of  fruit  of 
age  class  to  be  damaged  is 
greater  than  or  equal  to 
number  available,  then 
set  VDMG  = 1 .0  and  store 
number  of  any  remaining  fruit. 
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If  more  fruit  is  available  for 
this  age  class  than  needs  to 
be  damaged,  then  calculate 
percent  damaged  per  age 
class  (VDMG). 


For  each  plant  region  (IR), 
for  each  fruiting  site  on 
fruiting  branch  (M),  and 
for  each  fruit  age  class  (1C): 


Count  number  of  fruit 
damaged  (ACTDMG). 


If  number  of  fruit 
remaining  to  be  damaged 
after  considering  vegetative 
stems  is  greater  than  0 and 
number  of  available  fruit 
is  not  equal  to  remaining 
fruit  to  be  damaged,  then 
distribute  remaining  fruit 
to  be  damaged  evenly  over 
all  fruit  on  plant  (ADDMG). 
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This  next  section  calculates  actual 
percent  damage  for  each  fruiting  site. 


For  each  region  : 


For  each  sympodium  (L)  in  region  and 
each  fruiting  site  on  sympodium  (M)  : 
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For  each  fruit  age  class  : 


If  all  fruit  within  this  aqe  class 
is  to  be  damaged,  then  set 
percent  to  remove  (HZDMG) 
equal  to  1 .0. 


Skip  to  next  age  class. 


If  all  fruit  within  this  age  class 
in  this  region  is  to  be  damaged, 
then  set  percent  to  remove 
(HZDMG)  equal  to  1 .0. 
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If  all  fruit  in  this  region  of  this  age  class 
is  damaged,  set  H2DMG  = 1 . 

This  will  cause  all  fruit  within  this  region 
of  this  age  class  to  be  removed. 


Skip  to  next  age  class. 


Calculate  fraction  of  fruit  to  remove  from 
this  fruiting  site. 


If  fraction  to  remove  from  this  fruiting 
site  becomes  greater  than  1 , set  it  to  1 . 


Calculate  fraction  of  fruit  within  this 
age  class  to  remove  on  vegetative  stems. 
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Actual  removal  of  fruit  from  each 
fruiting  site  due  to  Heliothis  feeding. 


Is  fruit  a square? 


Square  weight  of  this  site  is 
decreased. 


Accumulate  amount  of  weight  loss. 


Reduce  fraction  of  fruit  remaininq 
at  this  site. 


If  fraction  of  fruit  at  site  goes 
to  less  than  0,  set  it  to  0. 
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Total  number  of  squares  on  plant 
is  decreased. 


Accumulate  amount  of  damage  on  each 
fruiting  site. 


Is  fruit  a green  boll? 


Calculate  weight  loss. 


Reduce  amount  of  seed  and  burr  nitrogen. 


Reduce  weight  of  green  boll. 


Decrease  total  number  of  green  bolls 
on  plant. 


Decrease  fraction  of  fruit  at  this  site. 


36 


FRUTLOS  = FRUIT  - FFRUT(K,L,M) 

BOL  ABZ  = BOL  ABZ  + FRUTLOS 

CHZDMG(K,L,M)  = CHZDMG(K,L,M)  + FRUTLOS 


If  fraction  of  fruit  at  site 
becomes  less  than  0,  set  it  to  0. 


Accumulate  damaged  fruit,  lost 
green  bolls,  and  amount  of  damage 
at  site. 


END 
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REDIST 


Calculate  fruit  remaining 
to  be  damaged  vithin  region. 


If  difference  is  less  than 
0.0001  or  no  fruit  is 
of  this  age  class  in  this  region  t 
skip  rest  of  logic. 


Subtract  number  to  be 
damaged  from  number 
already  marked  for  damage, 
and  calculate  remaining 
positions  in  region. 


If  remaining  positions 
become  0,  then  set  it  at 
1 .0. 


Calculate  proportion  to  be 
damaged  per  position. 


38 


A 


If  no  fruit  is  of  this  aqe 
class  at  this  position,  then 
set  N equal  to  last  position  with 
fruit  of  this  age  class. 


Calculate  number  of  fruit  to  be  damaged 
at  this  position. 


If  number  available  to  be  damaged  is  less 
than  number  to  be  damaged,  then  mark 
all  this  fruit  for  damage. 
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FRMTRX 


For  each  monopodium  : 


For  each  sympodium  : 


For  each  fruiting  site  on  sympodium  : 


Convert  proportion  of  fruit  damaged 
at  this  fruiting  site  to  a percent. 
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</></></></></>{/></>  4-1  </>  </>  </>  </>  0>  </>  </>  <0-  CO- 


PROGRAM  LISTING 


Command  Procedure  INMENU 

$ ON  CONTROL_Y  THEN  GOTO  END2 
$ OPEN/WRITE  OUTFILE  INPUT. DMG 
$ TYPE  SYS$ INPUT 

input  Heliothis  data 

1 . Heliothis  model 

2.  Scouting  report 

3 . No  damage 

$ INQUIRE  A "SELECT  NUMBER" 

$ IF  A .EQS.  1 THEN  GOTO  MODEL 

$ IF  A .EQS.  2 THEN  GOTO  SCOUT 

$ IF  A .EQS.  3 THEN  GOTO  NODMG 

$ MODEL: 

$ TYPE  SYS$ INPUT 


input  data  from  Heliothis  model 

A.  Enter  data  from  Heliothis  model 

B.  Enter  data  from  MOTHZV 

INQUIRE  A "SELECT  LETTER" 

IF  A .EQS.  "A"  THEN  GOTO  HELMODEL 
IF  A .EQS.  "B"  THEN  GOTO  MOTHZV 
HELMODEL: 

OPTION  = 1 

WRITE  OUTFILE  OPTION 
TYPE  SYS$ INPUT 

Enter  Heliothis  data  in  the  following  format 

(leave  1 space  between  numbers  - press  [RETURN]  on  blank  line  when 
inished) 

YEAR  DAY  NUMBER/ACRE  PHYSIOLOGICAL  AGE 

HEL_LOOP : 

INQUIRE  DATA 

IF  DATA  .EQS.  ""  THEN  GOTO  END 
WRITE  OUTFILE  DATA 
GOTO  HEL_LOOP 
MOTHZV: 

OPTION  = 2 

WRITE  OUTFILE  OPTION 
TYPE  SYS$ INPUT 

Enter  Hel iothis  data  in  the  following  format 

(leave  1 space  between  numbers  - press  [RETURN]  on  blank  line  when 
finished) 
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<j>  <r>  -co  <o  -co  <r>  <r>  </>  <r>  <n  Mi  <r>  <s>  <r>  <f>  <r>  </></>  <s> 


NUMBER  FRUIT  DAMAGED/ ACRE  BY 
YEAR  DAY  SMALL  LARVAE  LARGE  LARVAE 

MOTH_LOOP: 

INQUIRE  DATA 

IF  DATA  .EQS.  ""  THEN  GOTO  END 
WRITE  OUTFILE  DATA 
GOTO  MOTH_LOOP 
SCOUT: 

OPTION  = 3 

WRITE  OUTFILE  OPTION 
TYPE  SYS$ INPUT 


Enter  scouting  report  in  the  following  format 

(leave  1 space  between  numbers  - press  [RETURN]  on  blank  line  when 
inished) 


EAR  DAY  % SQUARE  DAMAGE  % BOLL  DAMAGE 

SCOUT_LOOP : 

INQUIRE  DATA 

IF  DATA  .EQS.  ""  THEN  GOTO  END 
WRITE  OUTFILE  DATA 
GOTO  SCOUT_LOOP 
NODMG: 

OPTION  = 4 

WRITE  OUTFILE  OPTION 

END: 

TYPE  SYS$ INPUT 

The  file  INPUT. DMG  has  been  created 

$ CLOSE  OUTFILE 
$ EXIT 
$ END2: 

$ CLOSE  OUTFILE 
$ EXIT 
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Subroutine  RDDMG 


SUBROUTINE  RDDMG 
£ ********************** 

C * READS  DAMAGE  FILE  * 

c 

INCLUDE  ' GOSCOM . FOR ' 

C 

LOGICAL  BUGGY 
C 

BUGGY  = .FALSE. 

HZBUG  = 0 
C 

REWIND  (50) 

C INPUT  HELIOTHIS  DATA 

READ (50,*)  I OPT I ON 
IF (IOPTION.EQ.1)  THEN 

10  READ (50  ,* ,END=15)  I BUGDAY , NUM , BUGAGE 
IF( (IBUGDAY-EMERGE) .EQ.IDAY)  THEN 

HZBUG  = NUM  * (0.035  + (0.093  * BUGAGE)) 

IF (HZBUG. LT.0)  HZBUG=0 
IF (BUGAGE. LE. 3)  NSTR=1 

IF(BUGAGE.GT.3.AND.BUGAGE.LE.4 .59)  NSTR=2 
IF (BUGAGE. GT. 4. 59. AND. BUGAGE. LE. 6. 79)  NSTR=3 
I F ( BUGAGE . GT . 6 . 79 . AND . BUGAGE . LE . 8 . 75 ) NSTR=4 
IF ( BUGAGE. GT. 8. 75)  NSTR=5 
BUGGY  = .TRUE. 

CALL  PREDMG 
END  IF 
GO  TO  10 
END  IF 

C ENTER  DATA  FROM  MOTHZV 

IF ( IOPTION.EQ. 2 ) THEN 

11  READ (50 ,* ,END=15)  I BUGDAY,  DMGSML , DMGLRGE 
IF ( (IBUGDAY-EMERGE) .EQ. IDAY)  THEN 

IF ( DMGSML. GT.0)  THEN 
NSTR  = 1 
HZBUG  = DMGSML 
CALL  PREDMG 
BUGGY  = .TRUE. 

END  IF 

IF (DMGLRGE. GT.0)  THEN 
NSTR  = 2 
HZBUG  = DMGLRGE 
CALL  PREDMG 
BUGGY  = .TRUE. 

END  IF 
END  IF 
GO  TO  11 
END  IF 
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15  IF (BUGGY)  CALL  DAMAGE 
C ENTER  SCOUTING  REPORT  DATA 
IF(I0PTI0N.EQ.3)  THEN 
25  READ (50 ,* ,END=30 ) I BUGDAY , SQP , GBP 
I F ( ( I BUGDAY-EMERGE ) . EQ . IDAY ) THEN 

I F ( SQP . GT . 0 . OR . GBP . GT . 0 ) CALL  PREDMG 
CALL  DAMAGE 
END  IF 
GO  TO  25 
END  IF 
30  RETURN 
END 
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Subroutine  PREDMG 


SUBROUTINE  PREDMG 

Q ******************************************************* 

C * DISTRIBUTES  NUMBER  OF  DAMAGED  FRUIT  ACROSS  AGE  * 

C * CLASSES  OF  FRUIT.  DISTRIBUTION  OF  DAMAGE  IS  BASED  * 

C * ON  INFORMATION  FROM  Wilson,  L.T.,  and  A.P.  Gutierrez* 

C * Fruit  predation  submodel:  Heliothis  larvae  feeding  * 

C * upon  cotton  fruiting  structures.  * 

C * Hilgardia  48(2):24-36,  1980  * 

Q ******************************************************* 

C 

INCLUDE  'GOSCOM.FOR ' 

C 

REAL  NUM 

DIMENSION  PROP (7) ,PDMG(5,7) ,PREF(5,7) ,STRAT(5,7) 

DIMENSION  PREFA(2,7) , STRATA(2,7),  TCLASS(7) 

DATA  PREF/. 6693,. 531,. 3855,. 105,. 0172, 1,1,. 8733,. 5641, 

& .3287,. 4124,. 4071,. 7412,. 8409,. 6701, 0,0, 1,1, 1,0,0, 

& .4119,. 3801,. 6018, 0,0,0,. 3783,. 6987, 0,0, 0,0,. 5498/ 

DATA  STRAT/1 , 1 , . 7014 , .479 , . 2068 , . 6485 , .9839 ,1,1,. 5064 , 

& .2876,  .4894, .6132, .9086, .5719,0,0, .2458, .9569,1,0,0, .0884 
& .2314,. 1944, 0,0,0,. 1004,. 14, 0,0, 0,0,. 031/ 

DATA  PREFA/. 5286,. 0611,. 9578,. 4464,. 5202,. 7555,. 3333, 1.0, 

& .1373, .491,0, .5385,0, .2749/ 

DATA  STRATA/. 9005,. 3429, .8775 , .7532 , .4634 , .7403 , .0819, .983 
& .0295, .2129,0, .1202,0, .0155/ 

SET  COUNTERS  TO  0 


20 

30 


SFRUIT=0 
TOTSQR  = 0 
TOTGB  = 0 
DENOM  = 0 
DO  30  1=1,7 

DO  20  J=1 ,5 
TCLASS ( I ) = 0 
PDMG (J, I ) = 0 
CONTINUE 
CONTINUE 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


IF  DAMAGE  IS  TO  OCCUR  TODAY,  CLASSIFY  FRUIT  AND  CALCULATE 
PROPORTION  OF  EACH  FRUIT  CLASS. 

CLASS  DESCRIPTION 

★★★★★★★★★★★★★★★★★★★★★★★★ 


SMALL  SQUARE 
MEDIUM  SQUARE 
LARGE  SQUARE 
FLOWER 

SMALL  GREEN  BOLL 
MEDIUM  GREEN  BOLL 
LARGE  GREEN  BOLL 


non 


DO  60  K=l,3 
DO  50  L=1 , 30 
DO  40  M=1 ,5 

CLASS (K,L,M)  = 0 
IF(FFRUT(K,L,M) .LE.0)  GO  TO  40 

IF(AGE(K,L,M) .LT. 2 . OR. AGE (K,L,M) .GT.48)  GO  TO  40 
IF(PCODE  (K, L,M)  .EQ.2.0R.FC0DE(K,L,M)  .GT.5)  THEN 
IF(AGEBOL (K,L,M) .LT.3)  CLASS (K,L,M) =4 
IF(AGEBOL(K,L,M) .GE. 3 . AND. AGEBOL (K,L,M) .LT.10) 

& CLASS (K,L,M) =5 

IF (AGEBOL (K,L,M) .GE. 10 .AND. AGEBOL (K,L,M) .LT.17) 
& CLASS (K , L , M)  =6 

IF  (AGEBOL  (K,L,M)  .GE.17) 

& CLASS (K,  L,M)  =7 

END  IF 

IF(FCODE(K,L,M) .EQ.l.OR.FCODE (K,L,M) .EQ.5)  THEN 
IF(AGE(K,L,M) .GE. 2 . AND. AGE (K, L,M) .LE.8) 

& CLASS  (K,L,M)=1 

IF(AGE(K,L,M) .GE. 9. AND. AGE (K,L,M) .LE.16) 

& CLASS (K, L,M) =2 

IF(AGE(K,L,M) .GE.17)  CLASS (K,L,M) =3 
END  IF 

IF (CLASS (K,L,M) .GT.0)  THEN 

SFRUIT  = SFRUIT  + FFRUT (K,L,M) 

IC  = CLASS (K,L,M) 

TCLASS(IC)  = TCLASS(IC)  + FFRUT  (K,L,M) 
IF(IC.LE.3)  TOTSQR  = TOTSQR  + FFRUT (K,L,M) 

IF ( IC.GE. 4)  TOTGB  = TOTGB  + FFRUT (K,Lf M) 

END  IF 
40  CONTINUE 

50  CONTINUE 
60  CONTINUE 

OPTION  = 1;  INPUT  NUMBER  OF  FRUIT  DAMAGED  BY  EACH  INSTAR 


*********************************************************************** 


* VARIABLES:  * 

* * 


* PREF( instar , fruit  class) 

* 

* PROP (fruit  class) 

* STRAT ( instar , fruit  class) 

* 

* PDMG ( instar , fruit  class) 

* 

* HZBUG 

* NSTR 


larval  food  preference  values  for  * 
various  age  fruits.  * 
proportion  of  each  age  class.  * 
age  specific  larval  stratification  * 
coefficients.  * 
proportion  of  damage  occurring  on  each  * 
fruit  age  class  by  each  instar.  * 
number  of  larvae/ acre.  * 
instar . * 
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* SMLARV (fruit  class)  = number  of  fruit  of  each  age  class/acre  * 

* damaged  by  larvae  from  instar  1-3.  * 

* LGLARV (fruit  class)  = number  of  fruit  of  each  age  class/acre  * 

* damaged  by  larvae  from  instar  4-5.  * 

*********************************************************************** 

C 

IF(lOPTION.EQ.l)  THEN 
DO  70  1=1,7 

IF(SFRUIT.LE.0)  RETURN 
PROP (I)  = TCLASS(I)  / SFRUIT 
70  CONTINUE 

DO  80  1=1,7 

80  DENOM  = PROP (I ) *PREF (NSTR , I ) *STRAT (NSTR, I ) + DENOM 

DO  90  1=1,7 

IF(DENOM.GT.0)  THEN 

PDMG (NSTR, I ) = (PROP (I ) *PREF (NSTR , I ) *STRAT (NSTR, I ) ) /DENOM 
END  IF 
90  CONTINUE 

DO  100  1=1,7 

IF (NSTR. LE. 3)  THEN 

SMLARV (I)  = HZBUG*PDMG (NSTR , I ) +SMLARV ( I ) 

ELSE 

LGLARV (I)  = HZBUG*PDMG (NSTR, I ) + LGLARV ( I ) 

END  IF 
100  CONTINUE 

END  IF 

I OPT I ON  = 2;  INPUT  DATA  FROM  MOTHZV 

IF(lOPTION.EQ.2)  THEN 
DO  110  1=1,7 

IF (SFRUIT. LE.0)  RETURN 
PROP (I)  = TCLASS ( I ) / SFRUIT 
110  CONTINUE 

DO  120  1=1,7 

120  DENOM  = PROP ( I ) *PREFA (NSTR , I ) *STRATA (NSTR , I ) +DENOM 
DO  130  1=1,7 

IF (DENOM. GT.0)  THEN 

PDMG (NSTR , I) = (PROP ( I ) *PREFA (NSTR , I ) *STRATA (NSTR , I ) ) /DENOM 
END  IF 
130  CONTINUE 

DO  140  1=1,7 

IF(NSTR.EQ.l)  THEN 

SMLARV (I)  = HZBUG*PDMG (NSTR , I ) +SMLARV ( I ) 

ELSE 

LGLARV (I)  = HZBUG*PDMG (NSTR , I ) +LGLARV ( I ) 

END  IF 
140  CONTINUE 

END  IF 
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I OPT I ON  = 3; 
ORE 


INPUT  PERCENT  OF  SQUARES  AND/OR  GREEN  BOLLS  DAMAGED  PER 


NUMBER  OF  FRUIT  OF  EACH  AGE  CLASS  TO  BE  DAMAGED  IS  DIRECTLY 
PROPORTIONAL  TO  NUMBER  OF  EACH  AGE  CLASS  PRESENT. 

IF(I0PTI0N.EQ.3)  THEN 
IF(SQP.GT.0)  THEN 
DO  150  1=1,3 

150  IF(TOTSQR.GT.0)  PROP(l)  = TCLASS(l)  / TOTSQR 

NUM  = TOTSQR* (SQP/100) *POPPLT 
DO  160  1=1,3 

160  FRUTDMG ( I ) = PROP(l)*NUM 

END  IF 

IF(GBP.GT.0)  THEN 
DO  170  1=4,7 

170  IF(TOTGB.GT.0)  PROP(l)  = TCLASS(l)  / TOTGB 

NUM  = TOTGB* (GBP/100) *POPPLT 
DO  180  1=4,7 

180  FRUTDMG (I)  = PROP(l)*NUM 

END  IF 
END  IF 
RETURN 
END 
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Subroutine  DAMAGE 


Q ********************************************************************* 

C * * 

C * DAMAGE  SUBROUTINE  * 

C * * 

C * THIS  SUBROUTINE  DISTRIBUTES  DAMAGE  CAUSED  BY  HELIOTHIS  * 

C * FEEDING.  * 

C * * 

(2  ********************************************************************* 
C 

INCLUDE  ' GOSCOM.FOR ' 

DIMENSION  BPER  (5)  , DTOTAL(7),  VEG(7),  BUGDMG  ( 7 ) , VDMG(7), 

& REGDMG  (4,7)  , BEGIN  (4)  , END  (4)  , DMG (7) 

INTEGER  UPPER 
C 

DATA  IPOS/0,0,0,5/ 

DATA  BPER/. 54,. 23,. 12,. 055,. 055/ 

SET  ALL  COUNTERS  TO  0 

SFRUIT  = 0 

ACTDMG  = 0 

CDMG  = 0 

REMN  = 0 

REMAIN  = 0 

ADDMG  = 0 

FRUTLOS  = 0 

HELIOTHIS  = 0 

PVEG  = 0 

VEGP  = 0 

C 

DO  30  1=1,7 
DO  20  K=l,4 
DO  10  L=1 ,5 


DTOTAL ( I ) 

= 0 

VEG(I) 

= 0 

BUGDMG (I) 

= 0 

VDMG(I) 

= 0 

REG (K,l) 

= 0 

CREG (K, I ) 

= 0 

REGDMG (K, I) 

= 0 

TCLS (K,L, I ) 

= 0 

POSDMG (K,L, I ) 

= 0 

10  CONTINUE 

20  CONTINUE 
30  CONTINUE 
C 

C DIVIDE  PLANT  INTO  THIRDS.  TOP  1/3  OF  PLANT  IS  REGION  1,  MIDDLE  1/3  IS 
C REGION J2,|AND  LOWER  1/3  IS  REGION  3.  REGION  4 IS  VEGETATIVE  STEM. 

C BOTH  VEGETATIVE  STEMS  ARE  GROUPED  UNDER  REGION  4. 
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c 

LOWER  = NFBR(l)  / 3.0 
MID  = LOWER  * 2.0 
UPPER  = NFBR(l) 

C 

C COUNT  TOTAL  NUMBER  OF  POSITIONS/ REG ION. 

C COUNT  NUMBER  OF  FRUIT/ AGE  CLASS  ON  MAIN  AND  VEGETATIVE  STEMS. 

C COUNT  NUMBER  OF  FRUIT/ AGE  CLASS/ REGION. 

C COUNT  NUMBER  OF  FRUIT/ AGE  CLASS/REGION/POSITION. 

C 

★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 

* VARIABLES:  * 

* * 


* SFRUIT 

* CLASS  (K,L,M) 

* 

★ 


* IPOS (IR) 

* DTOTAL (CLASS) 

* VEG (CLASS) 

* REG ( IR, CLASS) 

* TCLS (IR ,M, CLASS) 


# fruit  susceptible  for  damage. 

Age  class  of  fruit. 

1=1-8  4 = 25-26  7 = 41-48 

2 = 9-16  5 = 27-33 

3 = 17-24  6 = 34-40 

Max.  # of  posit ions/ region. 

# of  fruit/age  class  on  mainstem. 

# of  fruit/age  class  on  veg.  stems. 

# of  fruit/ age  class/region. 

# of  fruit/reg ion/posit  ion/age  class. 


* 

★ 

★ 

★ 

★ 

★ 

★ 

★ 

★ 

★ 


C 


DO  70  K=1 ,3 

IF(K.EQ.l)  THEN 
I START  = 1 
I END  = 3 

ELSE 

I START  = 4 
I END  = 4 

END  IF 

DO  60  I R= I START, I END 


IF(IR.EQ.I) 
BEGIN (IR) 
END ( IR) 
END  IF 
IF(IR.EQ.2) 
BEGIN (IR) 
END  (IR) 
END  IF 
IF(IR.EQ.3) 
BEGIN (IR) 
END (IR) 
END  IF 
IF(IR.EQ.4) 
BEGIN (IR) 
END (IR) 


THEN 

= MID  + 1 
= UPPER 

THEN 

= LOWER  + 
= MID 

THEN 
= 1 

= LOWER 

THEN 
= 1 
= 30 


1 
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END  IF 

DO  50  L=BEGIN(IR) ,END(IR) 

IF(BEGIN(IR) .EQ.0)  GO  TO  70 
DO  40  M=l,5 

IF (CLASS (K,L,M) .LE.0)  GO  TO  40 
ICLS  = CLASS (K,L,M) 

IF(K.EQ.l)  THEN 

DTOTAL (ICLS)  = DTOTAL (ICLS)  + FFRUT (K,L,M) 

IF(M.GT. IPOS ( IR) ) IPOS(IR)  = M 
ELSE 

VEG (ICLS)  = VEG (ICLS)  + FFRUT (K,L,M) 

IPOS ( IR)  = 5 
END  IF 

REG (IR, ICLS)  = REG ( IR , ICLS)  + FFRUT (K,L,M) 

TCLS(IR,M,ICLS)  = TCLS ( IR,M, ICLS)  + FFRUT (K,L,M) 

SFRUIT  = SFRUIT  + FFRUT (K,L,M) 

40  CONTINUE 

50  CONTINUE 

60  CONTINUE 
70  CONTINUE 
C 

C CONVERT  NUMBER  OF  FRUIT  TO  BE  DAMAGED  TO  PER  PLANT  BASIS. 

C WEIGHT  VERTICAL  DISTRIBUTION  OF  DAMAGE  BASED  ON  AGE  OF  LARVA. 

C Ramalho,  F.S.,  et  al.,  1984,  J.  Econ.  Ent.  77:591-594. 

C 

************************************************************************ 
* VARIABLES: 


* SMLARV (CLASS) 

* LGLARV (CLASS) 

* FRUTDMG (CLASS) 

* PLTDMG 

* CREG ( I R , CLASS ) 

* BUGDMG (CLASS) 

* PVEG 

* VEGP 


# of  fruit/ age  small  larvae  to  be  damaged  today. 

# of  fruit/age  large  larvae  to  be  damaged  today. 

# of  fruit/age  damaged  today. 

# of  fruit  to  be  damaged  per  plant. 

# fruit  to  be  damaged/ region/age  class. 

# of  fruit  to  be  damaged/age  class. 

proportion  of  fruit  of  each  age  class  on  veg.  stem, 
proportion  of  fruit  of  each  age  class  on  ma instem. 


k ★★★★★★★★★★★★★★★★★★★★★★ 


c 

DO  80  ICLS=1,7 

IF (DTOTAL ( ICLS) .GT.0.)  THEN 

PVEG  = VEG (ICLS)  / (DTOTAL ( ICLS)  + VEG (ICLS)) 

ELSE 

PVEG  = 0 
END  IF 

VEGP  = 1 - PVEG 
I F ( SMLARV ( ICLS) . GT . 0 ) THEN 

PLTDMG  = SMLARV (ICLS)  / POPPLT 

CREG (1, ICLS)  = CREG (1, ICLS)  + (PLTDMG  * 0.50)  * VEGP 
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80 


CREG (2 , ICLS)  = CREG (2 , ICLS)  + (PLTDMG  * 0.37) 

CREG (3 , ICLS)  = CREG (3, ICLS)  + (PLTDMG  * 0.13) 

CREG  (4, ICLS)  = PLTDMG  * PVEG 
END  IF 

IF(LGLARV(ICLS) .GT.0)  THEN 

PLTDMG  = LGLARV (ICLS)  / POPPLT 
CREG  (1 , ICLS)  = CREG (1, ICLS)  + (PLTDMG  * 0.39) 

CREG (2 , ICLS)  = CREG (2, ICLS)  + (PLTDMG  * 0.43) 

CREG (3 , ICLS)  = CREG (3 , ICLS)  + (PLTDMG  * 0.18) 

CREG (4/ ICLS)  = PLTDMG  * PVEG 
END  IF 

IF(FRUTDMG(ICLS) .GT.0)  THEN 

PLTDMG  = FRUTDMG ( ICLS)  / POPPLT 
CREG (1 , ICLS)  = CREG (1 , ICLS)  + (PLTDMG  * 0.45) 

CREG (2 , ICLS)  = CREG (2 , ICLS)  + (PLTDMG  * 0.40) 

CREG (3 , ICLS)  = CREG (3 , ICLS)  + (PLTDMG  * 0.15) 

CREG (4 , ICLS)  = PLTDMG  * PVEG 
END  IF 

BUGDMG (ICLS)  = BUGDMG ( I CLS ) + CREG (1 , ICLS) 
BUGDMG (ICLS)  = BUGDMG (ICLS)  + CREG (2, ICLS) 
BUGDMG (ICLS)  = BUGDMG (ICLS)  + CREG (3, ICLS) 
SMLARV ( ICLS)  = 0 
LGLARV (ICLS)  = 0 
FRUTDMG (ICLS)  = 0 
CONTINUE 


* VEGP 

* VEGP 


* VEGP 

* VEGP 

* VEGP 


* VEGP 

* VEGP 

* VEGP 


COMPARE  NUMBER  OF  EACH  AGE  CLASS  PRESENT  TO  NUMBER  OF 
EACH  AGE  CLASS  TO  BE  DAMAGED  TODAY. 

* VARIABLES: 

* 

* CDMG  = Cumulate  today's  damage. 

* POSDMG ( IR,M, IC)  = # fruit  to  be  damaged/region/pos it  ion/age  class. 

* BPER(M)  = % of  damage  occurr ing/position. 

* REMN  = # fruit  left  over  to  be  damaged  to  be  added 

* to  next  position. 

C 

DO  120  IC=1 ,7 
C 

C COMPARE  TOTAL  NUMBER  OF  FRUIT  TO  BE  DAMAGED/ AGE  CLASS  WITH  TOTAL 
C NUMBER  AVAILABLE.  IF  NUMBER  OF  FRUIT  AVAILABLE  OF  THAT  AGE  CLASS  IS 
C LESS  THAN  OR  EQUAL  TO  NUMBER  TO  BE  DAMAGED,  THEN  REMOVE  ALL  THIS 
C FRUIT. 

C ADD  REMAINDER  TO  VEGETATIVE  STEMS. 

C 

IF ( BUGDMG (IC) .EQ.0)  GO  TO  120 
IF(DT0TAL(IC) .LE. BUGDMG (IC) ) THEN 
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DMG(IC)  = 1.0 

CREG (4 , IC)  = BUGDMG(IC)  - DTOTAL(IC)  + CREG(4,IC) 

ELSE 

C 

C COMPARE  NUMBER  OF  FRUIT/AGE  CLASS/REGION  AVAILABLE  WITH  NUMBER  OF 
C FRUIT/ AGE  CLASS/REGION  TO  DAMAGE.  IF  LESS  FRUIT  IS  AVAILABLE 
C THAN  NEEDS  TO  BE  DAMAGED , THEN  REMOVE  ALL  FRUIT  OF  THAT  AGE  CLASS  FROM 
C REGION  AND  ADD  REMAINDER  TO  NEXT  REGION. 

C 

DO  110  IR=1 , 3 

IF(REG(IR,IC) .LE.0)  THEN 

CREG ( IR+1, IC)  = CREG ( IR+1 , IC)  + CREG(IR,IC) 

GO  TO  110 
END  IF 

IF(REG(IR,IC) .L2.CREG(IR, IC) ) THEN 
REGDMG ( IR, IC)  =1.0 

CREG (IR+1 , IC) =CREG (IR+1 , IC) + (CREG (IR, IC)  - REG(IR,IC)) 

CDMG  = REG ( IR, IC) 

GO  TO  110 
ELSE 
C 

C IF  NUMBER  OF  FRUIT  AVAILABLE  IS  GREATER  THAN  NUMBER  TO  BE  DAMAGED, 

C THEN  EVALUATE  EACH  FRUIT  ON  PLANT  AND  CALCULATE  AMOUNT  OF  FRUIT 
C TO  REMOVE  PER  FRUITING  SITE. 

C 

CDMG  = 0 
K=1 

DO  100  M=1 , IPOS  (IR) 

IF (BEGIN ( IR) .EQ.0)  GO  TO  110 
DO  90  L=BEGIN(IR) ,END(IR) 

C 

C WEIGHT  HORIZONTAL  DISTRIBUTION  OF  DAMAGE  BASED  ON  POSITION  ON  FRUITING 
C BRANCH. 

C 54%  POSITION  1 5.5%  POSITION  4 

C 23%  POSITION  2 5.5%  POSITION  5 

C 12%  POSITION  3 

C Ramalho,  F.S.,  et  al.,  1984,  J.  Econ.  Ent.  77:591-594 
C 

IF ( CLASS (K,L,M) .EQ. IC. AND.BUGDMG (IC) .GT.0)  THEN 
POSDMG ( IR,M, IC)  = CREG ( IR, IC)  * BPER(M) 

IF  NUMBER  OF  FRUIT  AVAILABLE  IS  LESS  THAN  OR  EQUAL  TO  NUMBER  TO  BE 
DAMAGED,  THEN  MOVE  REMAINDER  TO  BE  DAMAGED  TO  NEXT  POSITION  ON  THAT 
FRUITING  BRANCH. 

IF (TCLS (IR,M, IC) .LE .POSDMG (IR ,M, IC) ) THEN 
IF (TCLS ( IR,M+1 , IC) .GT.0)  THEN 

REMN  = POSDMG (IR,M,IC)  - TCLS (IR,M, IC) 

POSDMG (IR,M+1,IC)  = POSDMG (IR,M+1,IC)  + REMN 
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90 

100 


P0SDMG(IR,M, IC)  = TCLS(IR,M,IC) 
END  IF 
END  IF 
END  IF 
CONTINUE 

CDMG  = CDMG  + POSDMG ( IR , M, IC) 

CONTINUE 
END  IF 


IF  NOT  ENOUGH  FRUIT  HAS  BEEN  DAMAGED  WITHIN  THAT  REGION,  THEN 
REDISTRIBUTE  REMAINDER  EVENLY. 

IF(CDMG.LT.CREG(IR,IC) ) CALL  RED I ST 
110  CONTINUE 

END  IF 
120  CONTINUE 

DETERMINE  NUMBER  TO  BE  DAMAGED  ON  VEG  STEM. 

IF  AMOUNT  OF  FRUIT  AVAILABLE  TO  BE  DAMAGED  IS  LESS  THAN  NUMBER  TO  BE 
DAMAGED, THEN  REMOVE  ALL  FRUIT  OF  THAT  AGE  CLASS  AND  STORE  REMAINDER. 
IF  THERE  IS  ENOUGH  FRUIT  AVAILABLE,  THEN  CALCULATE  PERCENT  TO  BE 
REMOVED  FROM  EACH  FRUITING  SITE. 

DO  130  IC=1 ,7 

IF (VEG ( IC) .GT.0 .0 .OR.CREG (4 , IC) .GT.0.0)  THEN 
IF(CREG (4, IC) .GE.VEG (IC) ) THEN 
VDMG(IC)  =1.0 

REMAIN  = REMAIN  + (CREG(4,IC)  - VEG(IC)) 

ELSE 

VDMG(IC)  = CREG  (4 , IC)  / VEG(IC) 

END  IF 
ELSE 

REMAIN  = REMAIN  + CREG(4,IC) 

END  IF 
130  CONTINUE 

COUNT  ACTUAL  NUMBER  OF  FRUIT  DAMAGED  TODAY. 

ACTDMG  = # fruit  actually  damaged  today. 

DO  160  IR=1 ,4 

DO  150  M=1 , IPOS  ( IR) 

DO  140  101,7 

ACTDMG  = ACTDMG  + POSDMG (IR ,M, IC) 

140  CONTINUE 

150  CONTINUE 
160  CONTINUE 


IF  ANY  FRUIT  IS  LEFT  OVER  TO  BE  DAMAGED  AND  NUMBER  TO  BE  DAMAGED  IS 
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C NOT  EQUAL  TO  NUMBER  OF  SUSCEPTIBLE  FRUIT  ON  PLANT,  THEN 
C REDISTRIBUTE  REMAINDER  TO  REST  OF  FRUIT. 

C 

* ADDMG  = # fruit  to  be  damaged  to  add  to  each  fruiting  site. 
C 

IF ( REMAIN .GT . 0 . AND . SFRUI T . NE . ACTDMG)  THEN 
ADDMG  = REMAIN  / (SFRUIT  - ACTDMG) 

END  IF 
C 

DO  210  K=1 ,3 

IF(K.EQ.l)  THEN 
I START  = 1 
I END  = 3 
ELSE 

I START  = 4 
I END  = 4 


END  IF 

DO  200  IR= I START, I END 

IF(BEGIN(IR) .EQ.0)  GO  TO  210 
DO  190  L=BEGIN(IR) ,END(IR) 

DO  180  M=1 , IPOS ( IR) 

IF(I/RUT(K,L,M) .LE.0)  GO  TO  180 

DO  170  IC=1 ,7 

IF (CLASS (K,L,M) .NE.IC)  GO  TO  170 


CALCULATE  AMOUNT  OF  ABSCISSION. 

********************************************************************** 

* VARIABLES:  * 

* * 

* HZDMG  = % to  remove  from  site  due  to  Heliothis  damage.  * 

* VDMG(IC)  = # fruit  to  be  damaged/age  class  on  veg.  stems.  * 

********************************************************************** 

C 

IF(K.EQ.l)  THEN 

IF(DMG(IC) .EQ.1.0)  THEN 
HZDMG  =1.0 
GO  TO  175 
END  IF 

IF(REGDMG(IR,IC) .EQ.1.0)  THEN 
HZDMG  =1.0 
GO  TO  175 
END  IF 

IF (TCLS (IR,M, IC) .GT.0.0)  THEN 

HZDMG  = (POSDMG (IR,M, IC)  / TCLS (IR,M, IC) ) + ADDMG 
IF ( HZDMG. GT. 1.0)  HZDMG  = 1.0 
END  IF 
ELSE 
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HZDMG  = VDMG(IC)  + ADDMG 
IF(HZDMG.GT.1.0)  HZDMG  = 1.0 
END  IF 

C CALCULATE  AMOUNT  OF  FRUIT  DAMAGED  PER  SITE. 

175  CONTINUE 

IF(HZDMG.GT.0)  THEN 

IF(FCODE(K,L,M) .EQ. 1 .OR. PCODE (K,L,M) .EQ.5)  THEN 
WTLOS  = SQRWT (K,L,M)  * HZDMG 
SQRWT (K,L,M)  = SQRWT (K,L,M)  - WTLOS 
PQFLR  = PQFLR  + WTLOS 
FRUIT  = FFRUT(K,L,M) 

FFRUT (K,L,M)  = FFRUT (K,L,M)  - (FFRUT (K,L,M) *HZDMG) 
IF(FFRUT(K,L,M) .LT.0)  FFRUT (K,L,M)  =0 
FRUTLOS  = FRUIT  - FFRUT (K,L,M) 

SQRZ  = SQRZ  - FRUTLOS 

CHZDMG (K,Lf M)  = CHZDMG (K,L,M)  + FRUTLOS 
END  IF 

IF(FCODE(K,L,M) .EQ.2.0R.FC0DE (K,LfM) .GT.5)  THEN 
WTLOS  = BOLWGT(K,L,M)  * HZDMG 
SEEDN  = SEEDN  - (WTLOS/ (COTXX+GBOLWT) ) *SEEDN 
BURRN  = BURRN  - (WTLOS/ (COTXX+GBOLWT) ) *BURRN 
BOLWGT ( K , L , M)  = BOLWGT (K, L,M)  - WTLOS 
GBOLWT  = GBOLWT  - WTLOS 
FRUIT  = FFRUT (K,L,M) 

FFRUT (K,L,M)  = FFRUT (K,L,M)  - (FFRUT (KfL,M) *HZDMG) 
IF(FFRUT(KfLfM) .LT.0)  FFRUT (K,L,M)  =0 
FRUTLOS  = FRUIT  - FFRUT (K,L,M) 

BOLABZ  = BOLABZ  + FRUTLOS 
CHZDMG (K,L,M)  = CHZDMG (K,L,M)  + FRUTLOS 
END  IF 
END  IF 
HZDMG  = 0 

170  CONTINUE 

180  CONTINUE 

190  CONTINUE 

200  CONTINUE 
210  CONTINUE 
RETURN 
END 
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o n 


Subroutine  REDIST 


*************************************************************** 


* * 

* SUBROUTINE  REDIST  * 

* * 

* THIS  SUBROUTINE  REDISTRIBUTES  REMAINING  DAMAGE.  * 

* * 


*************************************************************** 

C 

INCLUDE  'GOSCOM.FOR ' 


N = 0 

5 CONTINUE 

XDIFF  = CREG(IR,IC)  - CDMG 
IF (XDIFF.LT. 0.0001. OR. REG (IR, IC) .EQ.0.0)  GO  TO  15 
RDMG  = CREG (IR, IC)  - CDMG 
RPOS  = IPOS ( IR)  - N 

IF (RPOS. EQ.0.0)  RPOS  = 1.0 
ADDMG  = RDMG  / RPOS 
C 

CDMG  =0.0 

DO  10  M=1 , IPOS ( IR) 

IF(TCLS(IR,M,IC) .LE.0.0)  THEN 
N = M 
GO  TO  10 
END  IF 

P0SDMG(IR,M, IC)  = POSDMG ( IR , M, IC)  + ADDMG 
IF(TCLS(IR,M,IC)  .LT. POSDMG (IR,M,IC) ) THEN 
IF(TCLS (IR,M, IC) .GT.0.0)  THEN 

REMN  = POSDMG (IR,M,IC)  - TCLS ( IR,M, IC) 

POSDMG (IR,M+1,IC)  = POSDMG (IR,M+1,IC)  + REMN 
POSDMG (IR,M,IC)  = TCLS ( IR,M, IC) 

REMN  = 0 
GO  TO  15 
END  IF 
END  IF 

CDMG  = CDMG  + POSDMG (IR,M, IC) 

10  CONTINUE 
GO  TO  5 
15  CDMG  = 0 
ADDMG  = 0 
N = 0 
RETURN 
END 


58 


Subroutine  FRMTRX 


C 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

104 

110 

202 

302 

402 

502 

C 

C 


C 


10 

20 

30 

C 


40 


★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★a- 
★ * 

* THIS  SUBROUTINE  DISPLAYS  OUTPUT  OF  INITIATION,  BLOOM,  AND  * 

* OPEN  BOLL  DATE,  PERCENT  FRUIT  AT  EACH  SITE,  AND  AMOUNT  OF  * 

* DAMAGE  AT  EACH  SITE.  THIS  IS  MODIFICATION  OF  * 

* SUBROUTINE  MATURE.  * 

* ★ 

★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★a-** 


INCLUDE  ' GOSCOM . FOR ' 
DIMENSION  CHZDMGP (3,30,5) 


rORMATS 

FORMAT (' 1' ,15X, 'Percentage  of  Fruit  Damaged  at  Each', 

.'  Location',/) 

FORMAT (15X, 'Ma instem  (Kl) ' ,15X, 'Vegetative  Branch  1 (K2 ) ' / 
. 10X , ' Ml  M2  M3  M4  M5',8X,'M1  M2  M3  M4 
FORMATC  L'  ,13,5  (2X,F4.0)  ,4X,5  (2X,F4.0) ) 

FORMATC  L' ,13,5  (2X,F4.0) ) 

FORMATC  L' ,13,5  (2X,F4.0)  ,7X, 'Vegetative  Branch  2 (K3)') 
FORMATC  L' ,13,5  (2X,F4.0)  ,7X,'M1  M2  M3  M4  M5  ' ) 


DO  30  K=1 ,NVBRCH 
NBRCH  = NFBR(K) 

DO  20  L=l, NBRCH 
NNID  = NNOD (K,L) 

DO  10  M=1,NNID 
PERCENT  DAMAGE 

CHZDMGP (K,L,M)  = CHZDMG (K,L,M)  * 100.0 
CONTINUE 
CONTINUE 
CONTINUE 

PRINT  PERCENT  DAMAGE  PER  SITE 
WRITE (LSTNG, 104) 

WRITE (LSTNG, 110) 

DO  40  L=1 ,15 
LL  = L-10 


IF(L.LT.6) 

IF(L.EQ.6) 
IF(L.EQ.7) 
IF (L. EQ. 8 ) 
IF(L.EQ.9) 
IF (L. EQ. 10 ) 
IF(L.GT.10) 


WRITE (LSTNG, 202 )L, (CHZDMGP (1 ,L,M)  ,M=1,5) 
, (CHZDMGP (2, L,M)  ,M=1,5) 
WRITE  (LSTNG, 302 )L, (CHZDMGP (1,6, M)  ,M=1,5) 
WRITE (LSTNG, 302 )L, (CHZDMGP (1,7, M)  ,M=1,5) 
WRITE (LSTNG, 402 )L, (CHZDMGP (1,8, M)  ,M=1,5) 
WRITE (LSTNG, 502) L , (CHZDMGP (1 ,9 ,M)  ,M=1,5) 
WRITE (LSTNG,302)L, (CHZDMGP (1 , 10 ,M) ,M=1,5) 
WRITE (LSTNG, 202) L, (CHZDMGP (1,L,M)  ,M=1,5) 
, (CHZDMGP ( 3, LL,M) ,M=1,5) 


M5 ' ) 


CONTINUE 

RETURN 

END 


*U.S.  GOVERNMENT  PRINTING  OFFICE:  1 988— 5 16-000/81086 
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